Search for: All records

We performed radio searches for the “spider” millisecond pulsar (MSP) candidates 4FGL J0935.3+0901, 4FGL J1627.7+3219, and 4FGL J2212.4+0708 using the Green Bank Telescope in an attempt to detect the proposed radio counterpart of the multi-wavelength variability seen in each system. We observed using the VEGAS spectrometer, centered predominantly at 2165 MHz; however, we were also granted observations at 820 MHz for 4FGL J1627.7+3219. We performed acceleration searches on each data set usingPRESTOas well as additional jerk searches of select observations. We see no evidence of a radio counterpart in any of the observations for each of the three systems at this time. Additional observations, perhaps at different orbital phases (e.g., inferior conjunction), may yield detections of an MSP in the future. Therefore, we urge continued monitoring of these systems to fully characterize the radio nature, however faint or variable, of each system.

 

We simulate scattering delays from the interstellar medium to examine the effectiveness of three estimators in recovering these delays in pulsar timing data. Two of these estimators use the more traditional process of fitting autocorrelation functions to pulsar dynamic spectra to extract scintillation bandwidths, while the third estimator uses the newer technique of cyclic spectroscopy on baseband pulsar data to recover the interstellar medium’s impulse response function. We find that either fitting a Lorentzian or Gaussian distribution to an autocorrelation function or recovering the impulse response function from the cyclic spectrum are, on average, accurate in recovering scattering delays, although autocorrelation function estimators have a large variance, even at high signal-to-noise ratio (S/N). We find that, given sufficient S/N, cyclic spectroscopy is more accurate than both Gaussian and Lorentzian fitting for recovering scattering delays at specific epochs, suggesting that cyclic spectroscopy is a superior method for scattering estimation in high-quality data.

 

We investigate the effectiveness of the statistical radio frequency interference (RFI) mitigation technique spectral kurtosis ($\widehat{\mathit{SK}}$) in the face of simulated realistic RFI signals.$\widehat{\mathit{SK}}$estimates the kurtosis of a collection ofMpower values in a single channel and provides a detection metric that is able to discern between human-made RFI and incoherent astronomical signals of interest. We test the ability of$\widehat{\mathit{SK}}$to flag signals with various representative modulation types, data rates, duty cycles, and carrier frequencies. We flag with various accumulation lengthsMand implement multiscale$\widehat{\mathit{SK}}$, which combines information from adjacent time-frequency bins to mitigate weaknesses in single-scale$\widehat{\mathit{SK}}$. We find that signals with significant sidelobe emission from high data rates are harder to flag, as well as signals with a 50% effective duty cycle and weak signal-to-noise ratios. Multiscale$\widehat{\mathit{SK}}$with at least one extra channel can detect both the center channel and sideband interference, flagging greater than 90% as long as the bin channel width is wider in frequency than the RFI.

 

The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency,ν) nor a change in dispersion measure alone (which would produce a delay proportional toν⁻²). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015.

 

We present an occultation study of compact radio sources by the plasma tail of interstellar comet 2I/Borisov (C/2019 Q4) both pre- and near-perihelion using the Arecibo and Green Bank radio telescopes. The interplanetary scintillation technique was used to probe the plasma tail at thePband (302–352 MHz), 820 MHz, and theLband (1120–1730 MHz). The presence and absence of scintillation at different perpendicular distances from the central axis of the plasma tail suggests a narrow tail of less than 6′ at a distance of ∼10′ (∼10⁶km) from the comet nucleus. Data recorded during the occultation of B1019+083 on 2019 October 31 with the Arecibo Telescope covered the width of the plasma tail from its outer region to the central axis. The systematic increase in scintillation during the occultation provides the plasma properties associated with the tail when the comet was at its pre-perihelion phase. The excess level ofL-band scintillation indicates a plasma density enhancement of ∼15–20 times that of the background solar wind. The evolving shape of the observed scintillation power spectra across the tail from its edge to the central axis suggests a density spectrum flatter than Kolmogorov and that the plasma density irregularity scales present in the tail range between 10 and 700 km. The discovery of a high-frequency spectral excess corresponding to irregularity scales much smaller than the Fresnel scale suggests the presence of small-scale density structures in the plasma tail, likely caused by interaction between the solar wind and the plasma environment formed by the comet.

 

Terzan 5 is a rich globular cluster within the galactic bulge containing 39 known millisecond pulsars, the largest known population of any globular cluster. These faint pulsars do not have sufficient signal-to-noise ratio (S/N) to measure reliable flux density or polarization information from individual observations in general. We combined over 5.2 days of archival data, at 1500 and 2000 MHz, taken with the Green Bank Telescope over the past 12 years. We created high-S/N profiles for 32 of the pulsars and determined precise rotation measures (RMs) for 28. We used the RMs, pulsar positions, and dispersion measures to map the projected parallel component of the Galactic magnetic field toward the cluster. The 〈B_∣∣〉 shows a rough gradient of ∼6 nG arcsec⁻¹(∼160 nG pc⁻¹) or, fractionally, a change of ∼20% in the R.A. direction across the cluster, implying Galactic magnetic field variability at sub-parsec scales. We also measured average flux densitiesS_νfor the pulsars, ranging from ∼10μJy to ∼2 mJy, and an average spectral indexα= −1.35, whereS_ν∝ν^α. This spectral index is flatter than most known pulsars, likely a selection effect due to the high frequencies used in pulsar searches to mitigate dispersion and scattering. We used flux densities from each observation to constrain the scintillation properties toward the cluster, finding strong refractive modulation on timescales of months. The inferred pulsar luminosity function is roughly power law, with slope$(d\log N)/(d\log L)=-1$at the high-luminosity end. At the low-luminosity end, there are incompleteness effects, implying that Terzan 5 contains many more pulsars.

 

Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.

 

Analyses of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nanohertz frequency band. The most plausible source of this background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for this background and assess its significance make several simplifying assumptions, namely (i) Gaussianity, (ii) isotropy, and most often, (iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated data sets. The data-set length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15 yr data set. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated data sets, even though their fundamental assumptions are not strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.

 

Abstract Radio pulsar signals are significantly perturbed by their propagation through the ionized interstellar medium. In addition to the frequency-dependent pulse times of arrival due to dispersion, pulse shapes are also distorted and shifted, having been scattered by the inhomogeneous interstellar plasma, affecting pulse arrival times. Understanding the degree to which scattering affects pulsar timing is important for gravitational-wave detection with pulsar timing arrays (PTAs), which depend on the reliability of pulsars as stable clocks with an uncertainty of ∼100 ns or less over ∼10 yr or more. Scattering can be described as a convolution of the intrinsic pulse shape with an impulse response function representing the effects of multipath propagation. In previous studies, the technique of cyclic spectroscopy has been applied to pulsar signals to deconvolve the effects of scattering from the original emitted signals, increasing the overall timing precision. We present an analysis of simulated data to test the quality of deconvolution using cyclic spectroscopy over a range of parameters characterizing interstellar scattering and pulsar signal-to-noise ratio (S/N). We show that cyclic spectroscopy is most effective for high S/N and/or highly scattered pulsars. We conclude that cyclic spectroscopy could play an important role in scattering correction to distant populations of highly scattered pulsars not currently included in PTAs. For future telescopes and for current instruments such as the Green Bank Telescope upgraded with the ultrawide bandwidth receiver, cyclic spectroscopy could potentially double the number of PTA-quality pulsars. 

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05 yr period and low redshift (∼0.02) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using pulsar timing array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it to simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5 yr data set employing PTA signal models containing Earth term-only as well as Earth+pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricitye₀< 0.5 and the symmetric mass ratioη> 0.1. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of 88.1 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches fore₀< 0.5 andη> 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 10⁹and (1.81 ± 0.01) × 10⁹M_☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.